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Differential cross sections for elastic scattering of photons from the deuteron have recently been 
measured at the Tagged-Photon Facility at the MAX IV Laboratory in Lund, Sweden. These first 
new measurements in more than a decade further constrain the isoscalar electromagnetic polariz¬ 
abilities of the nucleon and provide the first-ever results above 100 MeV, where the sensitivity to 
the polarizabilities is increased. We add 23 points between 70 and 112 MeV, at angles 60°, 120° and 
150°. Analysis of these data using a Chiral Effective Field Theory indicates that the cross sections 
are both self-consistent and consistent with previous measurements. Extracted values of a s = [12.1 
± 0.8(stat) ± 0.2(BSR) ± 0.8(th)[ x 10 -4 fm 3 and ft = [2.4 ± 0.8(stat) ± 0.2(BSR) ± 0.8(th)[ 
x 10 -4 fm 3 are obtained from a fit to these 23 new data points. This paper presents in detail the 
experimental conditions and the data analysis used to extract the cross sections. 

PACS numbers: 25.20.Dc, 24.70.+S 

Keywords: Compton Scattering; Deuterium; Polarizabilities. 


I. INTRODUCTION 

A new result on the extraction of the nucleon elec¬ 
tromagnetic polarizabilities was recently reported, based 
on recent measurements of Compton scattering from the 
deuteron ,l|. This paper presents in detail the motiva¬ 
tion, configuration, data analysis and critical evaluation 
of the experiment reported therein, as well as a parame¬ 
ter extraction using only the new data. 

Low-energy nuclear Compton scattering jX —> 7 X ex¬ 
plores how the internal degrees of freedom of the target 
behave in the electric and magnetic fields of a real exter¬ 
nal photon. Since these fields induce radiation multipoles 
by displacing the target constituents, the energy depen¬ 
dence of the emitted radiation provides a stringent test 
of the symmetries and strengths which govern the inter¬ 
actions of the constituents with each other and with the 
photon; see e.g. a recent review 
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The proton response can be measured directly and 
cleanly using a 4 H target. The neutron, however, is much 
more difficult to study because it is unstable outside the 
nucleus and its coupling to photons is much weaker. Em¬ 
bedding the neutron into a stable nucleus allows its two- 
photon response to be reconstructed. An added benefit 
of this approach is that the signal from the neutron is 
enhanced through its interference with the contributions 
from the charged proton. It also enables one to probe how 
the photons couple to the charged pion-exchange currents 
which provide the bulk of nuclear binding. The deuteron, 
which is the simplest stable few-nucleon system, is an 
ideal target for Compton-scattering experiments as it 
provides a conceptually clean probe of our understand¬ 
ing of both single-hadron and few-nucleon physics at low 
energies. 

After subtracting binding effects, theorists utilize such 
data to extract the two-photon response of the individual 
nucleon to the static fields; see [2, 0] for details. First, 
one subtracts the Powell amplitudes [4| for photon scat¬ 
tering on a point-like spin- 1/2 nucleon with an anoma¬ 
lous magnetic moment. The remainder is then expanded 
into energy-dependent radiation multipoles of the inci¬ 
dent and outgoing photon fields. Finally, these coeffi¬ 
cients are extrapolated to the values at zero photon fre¬ 
quency w. In that limit, the leading contributions are 
quadratic in to. Their coefficients are called the static 
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electric dipole polarizability a^i and the static magnetic 
dipole polarizability (3mi and can be separated by differ¬ 
ent angular dependences. They parametrize the stiffness 
of the nucleon against El —> El and Ml —»■ Ml transi¬ 
tions at zero photon energy, respectively. 

A host of information about the hadron response is 
thus compressed into cxei and (3 mi, often referred to as 
“the polarizabilities”. They are experimentally not di¬ 
rectly accessible since assumptions about the energy de¬ 
pendence and conventions on how to separate one- and 
two-photon physics enter. Nonetheless, they summarize 
information on the entire spectrum of nucleonic excita¬ 
tions. By comparing the quantities extracted from data 
with fully dynamical lattice QCD extractions which are 
anticipated in the near future HH3, the polarizabilities 
will offer a stringent test of our understanding of Quan¬ 
tum Chromodynamics (QCD). Most notable is the op¬ 
portunity to explore the two degrees of freedom with the 
lowest excitation energy, namely the pion cloud around 
the nucleon and the A(1232) excitation. Since both of 
these are dominated by isospin-symmetric interactions, 
differences between proton and neutron polarizabilities 
signal the breaking of isospin and chiral symmetry, in 
concert with such effects from short-distance physics. Be¬ 
sides being fundamental nucleon properties, a^i and [3 mi 
also play a role in theoretical studies of the Lamb shift 
of muonic hydrogen and of the proton-neutron mass dif¬ 
ference, and dominate the uncertainties of both mm* 

As recently reviewed in Ref. [U, a statistically con¬ 
sistent proton Compton-scattering database contains a 
cornucopia of points between 30 and 170 MeV, with 
good angular coverage and statistical uncertainties usu¬ 
ally around 5%. Based on this extensive set, McGovern 
et al. [l 2 j extracted the proton polarizabilities in Chiral 
Effective Field Theory (%EFT), the extension of Chiral 
Perturbation Theory to include baryons, as 1 

a p = 10.65 ± 0.35(stat) ± 0.2(BSR) ± 0.3(th) 

j3 p = 3.15 =F 0.35(stat) ± 0.2(BSR) ± 0.3(th), 

with x 2 = 113.2 for 135 degrees of freedom. Theoretical 
uncertainties are separated from those induced by appli¬ 
cation of the Baldin Sum Rule (BSR) for the proton jl3l ]. 
The BSR is a variant of the optical theorem which uses 
proton photoabsorption cross-section data to provide the 
constraint that a p + (3 P = 13.8 ± 0.4 fl4l ]. 

By contrast, the neutron polarizabilities, as extracted 
from deuteron Compton scattering, are poorly deter¬ 
mined. The calculations related to neutron polarizabili¬ 
ties appear to be theoretically well under control [|J, but 
the experimental deuteron data are of smaller quantity 
and poorer quality than those of the proton. Three ex¬ 
periments have thus far constituted the world data: the 


1 We use the canonical units of 10 4 fm 3 for the nucleon polariz¬ 

abilities throughout. 


pioneering effort of Lucas et al. at 49 and 69 MeV [l5| ; 
the follow-up measurement by Lundin et al. (l6| which 
covered similar energies and angles; and the extension to 
95 MeV by Hornidge et al. E3- This statistically consis¬ 
tent database contains only 29 data points at 4 energies 
between 49 and 95 MeV, with limited angular coverage, 
typical statistical uncertainties of more than 7%, and typ¬ 
ical systematic uncertainties in excess of 4%. From these 
data, the isoscalar (average) nucleon polarizabilities were 
extracted using the same X'EFT methodology as for the 
proton as 

a s = 10.9 ± 0.9(stat) ± 0.2(BSR) ± 0.8(th) 

/3 S = 3.6 =F 0.9(stat) ± 0.2(BSR) ± 0.8(th), 

with x 2 = 24.2 for 25 degrees of freedom [2]. The re¬ 
sult is again constrained by a BSR for the nucleon. The 
isoscalar value 

a s +/3 s = 14.5 ± 0.4 (3) 

is found by combining the proton BSR above with em¬ 
pirical partial-wave amplitudes for pion photoproduc¬ 
tion on the neutron which lead to the neutron sum rule 
& n + /3 n = 15.2 ±0.4 |T8j. The uncertainty in the neutron 
BSR is highly correlated with that for the proton. 

Combining isoscalar and proton polarizabilities, 
Eqs. dU) and ([5]), leads to neutron values 

a n = 11.1 ± 1.8(stat) ± 0.2(BSR) ± 0.8(th) 
f3 n = 4.1 =f 1.8(stat) ± 0.2(BSR) ± 0.8(th), 

which are clearly dominated by the statistical uncertain¬ 
ties, which in turn are much larger than for the proton; 
see Eq. Jl]). 

An alternative extraction of neutron polarizabili¬ 
ties from the 7 data points measured in quasi-elastic 
2 H( 7,7 'n)p above 200 MeV [l9[ is consistent with these 
numbers. Again using the neutron BSR constraint, one 
finds: 

a n = 12.5 ± 1.8(stat)+J;g(sys) ^ 1-1 (th) 

(3 n = 2.7 ± 1.8(stat)l?;®(sys) ± l.l(th), 

where the theory uncertainty may be underesti¬ 
mated [2Cj. No extractions from heavier nuclei exist; 
good Compton data is available on 6 Li [2l], [22[ , but com¬ 
parable data have not been published for any other few- 
nucleon systems. A third technique, extracting the neu¬ 
tron polarizabilities from scattering from the Coulomb 
field of heavy nuclei, appears to be plagued by poorly 
understood systematic effects; see e.g. [2j. 

New deuteron data of good quality and reproducible 
systematic uncertainties are therefore necessary to see 
commonalities and differences in the two-photon re¬ 
sponses of protons and neutrons. The experiment de¬ 
tailed in this paper effectively doubled the deuteron 
Compton database and significantly reduced the uncer¬ 
tainties of the neutron polarizabilities. These data over¬ 
lap the previous sets at lower energies and add points up 
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TABLE I: Basic parameters of the electron beam and the 
tagged-photon beam for the two run periods, RP1 and RP2. 



RP1 

RP2 

-^'electron [MeV] 

144 

165 

E-f [MeV] 

65 - 97 

81 - 115 

£ 7 ,bin [MeV] 

69.6, 77.8 

85.8, 94.8 


86.1, 93.7 

103.8, 112.1 

AE 7! bin [MeV] 

~8.0 

~8.5 


to 112 MeV, with statistical and systematic uncertainties 
on par with preceding measurements. The extension to 
higher energies is particularly important since the sensi¬ 
tivity of the cross sections to the polarizabilities increases 
roughly with the square of the photon energy. 

The resulting augmented world database is statistically 
consistent and recently resulted in a new extraction of the 
isoscalar polarizabilities in Refs. eh as 

a s = 11.1 ± 0.6(stat) ± 0.2(BSR) ± 0.8(th) 

P s = 3.4 =f 0.6(stat) ± 0.2(BSR) q= 0.8(th). 

Our new measurement thus reduces the statistical uncer¬ 
tainties in a s and /3 S by a factor of 1/3. For the very first 
time, the uncertainty is now dominated by the theoretical 
uncertainties of the extraction. 

While some aspects of our findings have been summa¬ 
rized briefly in a recent publication (l||, we now provide 
a more detailed description of the entire experimental ef¬ 
fort, including complementary information to aid in the 
interpretation of the results. In Sections |TT| and IHII we 
present the experimental setup at MAX IV and the data- 
analysis process, and pay special attention to yield cor¬ 
rections and systematic uncertainties. Section IIVI con¬ 
tains the resulting cross sections and an extraction of the 
polarizabilities, focusing on the self-consistency of this 
data set and its agreement with previous measurements. 



FIG. 1: The layout of the experimental area showing the 
location of the tagging spectrometer, focal-plane liodoscope, 
deuterium target, and Nal(Tl) detectors labeled DIANA, 
BUNI, and CATS. 



FIG. 2: Enlarged diagram of the tagger magnet and FP 
hodoscope portion of the experimental layout (not to scale). 


II. EXPERIMENT 

The experiment was performed at the Tagged-Photon 
Facility mi located at the MAX IV Laboratory ]25j in 
Lund, Sweden. A pulse-stretched electron beam [26| with 
a typical current of 15 nA and a duty factor of 45% 
was used to produce quasi-monoenergetic photons via 
the bremsstrahlung-tagging technique |23, H • The basic 
parameters of the electron and resulting tagged-photon 
beams are given in Table Q] for the first and second run 
periods, RP1 and RP2. An overview of the experimental 
layout is shown in Fig. [l] 

The tagging magnet and focal-plane (FP) hodoscope 
[20t were used extensively at the Saskatchewan Accelera¬ 
tor Laboratory prior to their use at the MAX IV Labora¬ 
tory. The dipole field of the magnet is used to momentum 
analyze the post-bremsstrahlung electrons, which are de¬ 
tected in the FP hodoscope by 63 plastic scintillators. 


The scintillators are 25 mm wide and 3 mm thick and 
arranged into two rows. The rows are offset by 50% of 
the scintillator width, with each overlap defining a FP 
channel (see Fig. 0. The typical width of a FP chan¬ 
nel was ~400 keV and the nominal electron rate was 
1 MHz/channel. As the Compton counting rate is low, 
the focal plane was subdivided into four bins. The cen¬ 
tral photon energies for each bin, as well as the average 
bin width, are given in Table [I] 

The size of the photon beam was defined by a tapered 
tungsten-alloy primary collimator of 19 mm nominal di¬ 
ameter. The primary collimator was followed by a dipole 
magnet and a post-collimator which were used to remove 
any charged particles produced in the primary collimator. 
The beam spot at the target location was approximately 
60 mm in diameter. 

The tagging efficiency is the ratio of the number of 
tagged photons which struck the target to the number of 
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post-bremsstrahlung electrons which were registered by 
the associated FP channel. It depends on the collimator 
size and the electron-beam energy. It was measured dur¬ 
ing the experiment start-up with the three large-volume 
Nal(Tl) photon spectrometers placed directly in the low- 
intensity photon beam and was monitored during data 
collection on a daily basis via dedicated measurements 
with a compact lead-glass photon detector, which was 
easily raised into and lowered out of the photon beam. 

The liquid deuterium target used in this experi¬ 
ment was based on a design used in previous measure¬ 
ments 0 , but with measures developed to eliminate ice 
build-up on the target endcaps. These measures included 
a several-day bake-out of the vacuum vessel to reduce in¬ 
ternal gases, thicker Kapton foils for the vacuum cham¬ 
ber windows, and the implementation of a N 2 gas shield. 
These last two measures reduce the penetration of water 
vapor from the air into the insulation vacuum surround¬ 
ing the cell. The cell was cylindrical, 150 mm long and 
68 mm in diameter. The spherical endcaps were convex 
so that the total length of the cell was 170 mm. The cell 
was oriented so that its central axis was collinear with 
the beam line. The housing chamber was constructed of 
stainless steel with a thickness of ~1 mm in the vicinity 
of the scattering plane and ~2 mm elsewhere. 

Three large-volume, segmented Nal(Tl) detectors la¬ 
beled BUNI 0, CATS (HU, and DIANA H in Fig. □ 
were used to detect the Compton-scattered photons. The 
detectors were located at laboratory angles of 60°, 120°, 
and 150°. These detectors were each composed of a 
Nal(Tl) core surrounded by optically isolated, annular 
Nal(Tl) segments. The cores of the BUNI and CATS 
detectors each measure 26.7 cm in diameter, while the 
core of the DIANA detector measures 48.0 cm. The 
depth of all three detectors is greater than 20 radiation 
lengths. The annular segments are 11 cm thick on the 
BUNI and CATS detectors and 4 cm thick on the DI¬ 
ANA detector. Additionally, BUNI and CATS each has 
a plastic-scintillator annulus that surrounds the Nal(Tl) 
annulus. Each detector was shielded by lead with a front 
aperture that defined the detector acceptance. A plastic- 
scintillator paddle was placed in front of the aperture to 
identify and veto charged particles. The detectors have 
an energy resolution of better than 2% at energies near 
100 MeV. Such resolution is necessary to separate un¬ 
ambiguously the elastically scattered photons from those 
originating from the breakup of deuterium. 

The signals from each detector were passed to analog- 
to-digital converters (ADCs) and time-to-digital convert¬ 
ers (TDCs) and the data were recorded on an event-by¬ 
event basis. The experimental data were collected in two 
separate four-week run periods. The first run period em¬ 
ployed an electron-beam energy of 144 MeV; the second 
run period used 165 MeV. The first week of each period 
was dedicated to in-beam studies (see below), measure¬ 
ments of 12 C( 7 , 7 ) [13 to establish the absolute normal¬ 
ization and systematics of the setup, and cooling of the 
liquid-deuterium target. The remaining three weeks were 


used to perform the measurements on deuterium reported 
in this article. 


III. DATA ANALYSIS 

The Compton scattering cross section can be written 
as 

da Y , i i . . 

777 = (o- Jt;- JRJT, (7) 

ail Deff A 7 K e ff 

where (Y/f2 e ff) is the scattered-photon yield normalized 
to the effective solid-angle acceptance of the detector, N 7 
is the number of tagged photons incident upon the tar¬ 
get, and K e ff is the effective thickness of the target (the 
number of nuclei per unit area). /r and /t represent cor¬ 
rection factors due to rate- and target-dependent effects, 
respectively, as explained below. 

A. Acceptance-normalized yield 

1. Scattered-photon yield 

During the experiment, the ADCs allowed reconstruc¬ 
tion of the scattered-photon energies, while the TDCs 
enabled coincident timing between the Nal(Tl) detectors 
and the FP hodoscope. The energy calibration of each 
of the Nal(Tl) detectors was determined by placing it 
directly into the reduced-intensity photon beam and ob¬ 
serving its response as a function of tagged-photon en¬ 
ergy. To calibrate each ADC, a spectrum was filled (see 
the inset to Fig. m for each FP channel by selecting only 
on tagged photons for that channel. The position of the 
peak (in ADC channels) was then plotted against the ex¬ 
pected photon energy, as determined from the tagging 
magnet field map, for all the FP channels. As an exam¬ 
ple, the calibration for BUNI is shown in Fig. [3] 

Missing energy (ME) was defined to be the difference 
between the expected energy of the detected photon (as 
determined from the tagger magnet and FP hodoscope 
placement) and the energy deposited by the photon in the 
Nal(Tl) detector. In both BUNI and CATS, the energy 
deposition in the annulus was added to the core energy 
deposition to improve the resolution. The measured in¬ 
beam response for the BUNI detector, together with a 
fitted Geant4 (34J simulation of this response, is shown 
in Fig. [H The Geant4 simulation output was determined 
for the case of the BUNI detector positioned directly in 
the photon beam. This “intrinsic” simulation was then 
smeared with a Gaussian function according to 

( e ,- e 3 -- P2 ) 2 

/ t> inbeam _ \ ^ 1 „ -2 (Eavi\ z cinbeam /o\ 

Ki ~^ Ej p 3 ' ’ 1 

where E^ was the central energy of bin i(j ), »Sj nbeam was 
the number of counts in bin j of the simulated detector- 
response spectrum, and pi, 2,3 were fitting parameters 
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FIG. 3: A plot of the predicted photon energy for each 
FP channel versus the peak position of the ADC spectrum 
(Padc) for the same channel. The plot is fit with a linear 
function to determine the calibration of the Nal(Tl) detec¬ 
tor. (Inset) A typical tagged-photon ADC spectrum for BUNI 
showing the location of the ADC channel corresponding to the 
peak (Padc). 



FIG. 4: (Color online) In-beam detector response as a func¬ 
tion of missing energy for the BUNI detector. The data points 
are the result of summing the energy deposition in both the 
Nal(Tl) core and annulus. The solid red curve is the simulated 
Geant 4 detector response fitted to the data. The dashed line 
is the ME spectrum obtained by analyzing only the core crys¬ 
tal. The addition of the annulus energy improves the FWHM 
by almost 50%. 


that accounted for the individual characteristics of each 
Nal(Tl) detector, such as non-uniform doping of the crys¬ 
tal, that are difficult to model in Geant4. This process 
was repeated for CATS and DIANA. 

It was observed that the gain of the Nal(Tl) detectors 
drifted over the course of the run periods. In order to cor¬ 
rect for this gain drift, cosmic-ray data were collected im¬ 
mediately after each in-beam calibration run and prior to 


moving the detector to its scattering location. “Straight- 
through” cosmic rays, defined by requiring a large energy 
deposition in annular segments on opposite sides of the 
core crystal, were selected because these events should 
have a constant energy deposition in the detector. The 
gain drift for each PMT for each run could be determined 
by monitoring the location and shape of the resulting 
cosmic-ray peaks. These gain-drift corrections were then 
applied to the data. 

The energy calibration of the tagger focal plane was 
confirmed by observing highly energetic capture photons 
from the ir~ + d —> 7 + 2n reaction. This reaction chan¬ 
nel was present as the untagged bremssstrahlung spec¬ 
trum extended beyond pion production threshold energy, 
and the most probable energy of the capture photon is 
~131 MeV [35[. The agreement between the absolute 
photon energy and that reconstructed from the tagger 
FP energy calibration was better than 1%. 

Large backgrounds arose when the beam intensity 
was increased from 10-100 Hz per FP channel (for in¬ 
beam runs) to ~1 MHz per channel (for scattering 
runs). The dominant sources of background were un¬ 
tagged bremsstrahlung photons (which scaled with the 
beam intensity) and cosmic rays. These backgrounds 
obscured the timing and ME peaks of the elastic pho¬ 
tons in the TDC and ADC spectra, respectively. Cosmic 
rays deposit a large amount of energy in the detector 
overall and in the annular segments in particular. A cut 
placed on the Nal(Tl) annulus (BUNI and DIANA) or the 
plastic-scintillator (CATS) annulus removed ~95% of the 
cosmic-ray background from the scattering data. An ad¬ 
ditional cut placed on the thin plastic-scintillator paddle 
in front of each of the Nal(Tl) detectors removed charged 
particles. These cuts reduced the number of events by 
~50%. 

In order to further reduce the untagged 
bremsstrahlung background in the FP TDC spec¬ 
trum, a cut was placed on the energy deposited in 
the Nal(Tl) detector. Selecting only events with an 
energy deposition E min < E < E max , where U min ( max) 
is the minimum (maximum) tagged-photon energy, 
enabled the prompt peak to be identified in the FP TDC 
spectrum 2 (see Fig. [5]). This prompt peak represented 
coincidences between post-bremsstrahlung electrons in 
the FP hodoscope and elastically scattered photons in 
the Nal(Tl) detectors. 

For each Nal(Tl) detector and for each FP channel, 
events occurring within the prompt peak were selected 
and a prompt ME spectrum was filled. The process was 
repeated for a second cut placed on a purely accidental 
timing region, and an accidental spectrum was filled. 

For each run period, data collected during dedicated, 


2 The structure seen in the TDC spectrum was a result of the 
incomplete filling of the pulse-stretcher ring (T=108 ns) and the 
3.3 MHz extraction shaker (T=305 ns) [3611 . 
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FIG. 5: (Color online) The FP TDC spectrum for the scat¬ 
tering data. The prompt (red) and the accidental (shaded) 
windows are indicated. 


beam-off cosmic-ray runs were utilized to remove the 
~5% of cosmic-ray events that survived the annulus cuts. 
For each Nal(Tl) detector and for each FP channel, these 
data were subjected to the same annulus cuts above, 
and cosmic-ray spectra were filled. These spectra were 
scaled by the ratio of events with an energy exceeding the 
electron-beam energy in the prompt spectra to those in 
the cosmic-ray spectra, and then subsequently subtracted 
from the prompt spectra. The procedure was repeated 
for the accidental spectra. In this way, prompt and acci¬ 
dental spectra free from cosmic rays were produced. 

The Geant4 in-beam simulation was extended to reflect 
the experimental setup for scattering runs. A numerical 
function J- t was defined by 

j)=p 0 A+^r tter , (9) 

where Ai was the number of counts in bin i of the acciden¬ 
tal spectrum, po was the scale factor of the accidentals, 
and 7£| catter was given by Eq. ([8j using the scattering 
response spectrum from the Geant4 simulation (see the 
top panel of Fig. ©■ The range of the fitting window 
varied from as small as [-10,+10] MeV to as large as 
[ -20,+20] MeV in ME. By fitting several windows over 
this range, an estimate of the dependence of the extracted 
yield on the width of the fitting window was obtained. 
This kinematic-dependent uncertainty depended strongly 
on the signal-to-noise ratio in the prompt ME spectrum 
and ranged from 2% to 11%. 

The bottom panel of Fig.|B]shows a typical “true” scat¬ 
tering spectrum (prompts minus accidentals) together 
with the corresponding fitted, convoluted response func¬ 
tion ^scatter _ q^g sca ttering yield was extracted by in¬ 
tegrating the true spectrum over the region of interest 
(ROI) (—2.0 MeV < ME < 2.0 MeV) indicated by the 
vertical dashed lines. The ROI was not allowed to extend 
below —2.0 MeV as photons from the photodissociation 
of the deuteron are kinematically allowed in this region. 



FIG. 6: (Color online) Top panel: The result of fitting the 
sum (solid) of the accidental spectrum (dashed) and the fitted 
detector response to the prompt spectrum. Bottom panel: 
The true elastic peak (prompts minus accidentals) together 
with the fitted detector response (solid). The vertical, dashed 
lines indicate the region of integration over which the yield 
extraction occurs, —2.0 MeV < ME < 2.0 MeV. 


The elastically scattered photon yields determined ac¬ 
cording to this procedure are given in Table [TT] 

The —2.0 MeV < ME < 2.0 MeV ROI was care¬ 
fully chosen, taking into account the detector resolution 
(FWHM <2% at 100 MeV), to minimize the contribu¬ 
tion of photons from e?(7,7 , ) ri P to the extracted yield. 
Contamination from non-elastic photons was investigated 
with a direct and an indirect search. The direct method 
involved simulating the photons from the breakup reac¬ 
tion and adding this new lineshape to the fitting algo¬ 
rithm. This new lineshape was displaced 2.2 MeV to 
the low-energy side of the elastic peak to account for the 
reaction threshold. The data were re-fit and the contri¬ 
bution of the non-elastic photons to the extracted yield 
was calculated. This contribution was found to be consis¬ 
tent with zero within uncertainties. The indirect method 
relied on the quality of the fit shown in Fig. [Gj which has 
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a typical reduced y 2 value of < 1. If there were a sizable 
contribution arising from a non-elastic reaction, it would 
cause the extracted cross sections to vary depending on 
the ROI width. Since non-elastic photons lie to the left 
of the elastic peak, the lower edge of the ROI was varied 
by ±400 keV. It was found that the extracted cross sec¬ 
tions with the varied ROIs agreed with the ones listed in 
this paper within uncertainties. Thus, it was concluded 
that the cross sections presented here do not suffer from 
contributions from d^^^np. 


2. Detector Acceptance 

The detector acceptances depended on the cuts em¬ 
ployed during the data analysis and the width and loca¬ 
tion of the integration ROI. The geometrical solid angle 
subtended by the Nal(Tl) detectors was corrected for the 
geometry of the experimental setup. Both of these effects 
were studied using the Geant4 Monte Carlo simulations. 
The effective solid angle ft e s was given by 


n eS = 4tt, (io) 

-^events 

where N evenis was the total number of simulated events, 
and N ™ L was the number of simulated events that 
eventually populated the ROI. 

The cosmic-rejection cuts removed a very small num¬ 
ber of good Compton-scattered photons from the ROI. 
For the CATS (60°) detector, the plastic-scintillator an¬ 
nulus was used as a cosmic-ray veto. For the DIANA 
(150°) detector, the thin Nal(Tl) annulus was used to 
reject cosmic rays. In each case, the cosmic-ray veto was 
more than 20 cm from the central cylindrical symmetry 
axis of the detector and only ~1% of all elastically scat¬ 
tered photons were rejected. For the BUNI (120°) de¬ 
tector, annular Nal(Tl) segments ~13 cm from the cylin¬ 
drical symmetry axis of the detector were used for the 
rejection of cosmic rays. As a result, 6% of all scattered 
photons were rejected. However, the Nal(Tl) segments in 
BUNI provided sufficient energy resolution to determine 
the cut placement accurately resulting in a systematic 
uncertainty of ~2%. The charged-particle veto removed 
~1% of all scattered photons. 

The effective solid angle for each data point is given 
in Table [II] A sufficient number of events were simulated 
so that the statistical uncertainties are <1%. The sys¬ 
tematic uncertainties include the effects of the cosmic- 
rejection cuts as well as uncertainty from measurements 
of the target-detector distance and detector-aperture di¬ 
ameter, typically ±2 mm. 



FIG. 7: Time evolution of the tagging efficiency measured 
using the compact lead-glass detector for the E-, = 93.4 MeV 
bin from RP1. The shaded region indicates the uncertainty 
due to daily variations in £tag- 


B. Scale normalization 

1. Number of beam photons 

The number of beam photons incident on the target 
iV 7 was determined from 

Nj= N e - £ t ag, (11) 

where N e is the number of post-bremsstrahlung electrons 
detected in each FP channel and £ tag is the tagging effi¬ 
ciency. 

The number of electrons striking the focal plane (N e ) 
was counted by the FP scalers (see Table Hill) . The back¬ 
ground rate in the focal plane, obtained from beam-off 
runs, was on the order of 1 Hz per channel and was 
thus negligible compared to the beam-on electron rate of 
~1 MHz per channel. The tagging efficiency £ tag was de¬ 
termined from the ratio of the number of photons tagged 
by a FP channel and recorded in the in-beam photon 
detector to the number of post-bremsstrahlung electrons 
recorded by the same channel. Livetime-corrected beam- 
on and beam-off backgrounds (~5%) were removed from 
the data. The focal plane was divided into four bins, 
each 16 channels wide. The tagging efficiency for each 
bin was taken to be the electron-weighted average of the 
tagging efficiencies for each of the 16 channels. A plot of 
the daily tagging efficiency measured using the compact 
lead-glass detector for the E 1 = 93.4 MeV bin is shown 
in Fig. 0 

In order to observe any systematic difference between 
the tagging efficiency determined by the Nal(Tl) de¬ 
tectors and the lead-glass detector, measurements were 
taken with the CATS detector and immediately there¬ 
after with the lead-glass detector. These data allowed 
for a <2% correction to be made to the compact lead- 
glass detector results, mainly due to the larger volume 
of the Nal(Tl) and the presence of the thin paddle used 
to veto charged particles. This correction was applied 
to the average value of the tagging efficiency determined 









TABLE II: Extracted yields and effective solid angle at each energy and angle. For the yields, the first uncertainty is statistical 
and the second is a kinematic-dependent systematic. For the effective solid angles, the first uncertainty is statistical and the 
second is an angle-dependent systematic. The upper four energy bins are from RP1 and the lower four from RP2 (also in 
Tables ITTT1 ITVl andlVTl). 


(MeV) 

60° 

Y fleff (msr) 

120° 

Y D e ff (msr) 

150° 

Y Heff (msr) 

69.6 

77.9 

86.1 

93.4 

1080 ± 182 ± 61 29.0 ± 0.3 ± 1.2 

995 ± 137 ± 59 25.3 ± 0.3 ± 1.1 

809 ± 93 ± 18 28.1 ± 0.3 ± 1.2 

440 ± 63 ± 13 26.5 ± 0.3 ± 1.1 

1106 ± 191 ± 42 42.9 ± 0.4 ± 1.9 

1528 ± 152 ± 32 42.9 ± 0.4 ± 1.9 

1312 ± 115 ± 26 38.1 ± 0.4 ± 1.6 

1148 ± 90 ± 23 38.2 ± 0.4 ± 1.6 

1034 ± 142 ± 84 22.8 ± 0.2 ± 0.7 

790 ± 115 ± 41 22.6 ± 0.2 ± 0.7 

573 ± 90 ± 63 22.2 ± 0.2 ± 0.7 

85.8 

94.8 
103.8 
112.1 

1669 ± 199 ± 187 24.0 ± 0.2 ± 1.0 

1639 ± 161 ± 127 18.9 ± 0.2 ± 0.8 

1266 ± 117 ± 30 21.0 ± 0.2 ± 0.9 

842 ± 91 ± 21 21.5 ± 0.2 ± 0.9 

2199 ± 156 ± 187 41.6 ± 0.4 ± 1.8 

2633 ± 142 ± 94 41.1 ± 0.4 ± 1.8 

1919 ± 117 ± 55 39.2 ± 0.4 ± 1.7 

1370 ± 95 ± 29 37.7 ± 0.4 ± 1.6 

1616 ± 198 ± 50 20.2 ± 0.2 ± 0.6 

1587 ± 174 ± 48 19.8 ± 0.2 ± 0.6 

1424 ± 141 ± 71 19.3 ± 0.2 ± 0.6 

1034 ± 115 ± 29 19.4 ± 0.2 ± 0.6 


with the compact lead-glass detector. One-half of the 
correction was assigned as a systematic uncertainty. 

The number of post-bremsstrahlung electrons N e , tag¬ 
ging efficiencies £tag> and number of beam photons _/V 7 
are presented in Table uni 


2. Effective target thickness 


The effective target thickness n e s (the number of nuclei 
per unit area) was given by 


^eff — 


P L N A 
A 


( 12 ) 


where p was the average density of liquid deuterium, L 
was the effective target length, Na was Avogadro’s num¬ 
ber, and A was the molar mass (4.0282 g/mol). 

The target pressure and temperature were systemat¬ 
ically recorded for each run. As the density of liquid 
deuterium is related to its pressure [13] , an average tar¬ 
get density was calculated by determining the deuterium 
density measured during each run and weighting this den¬ 
sity by the number of FP electrons recorded in the same 
run. For each run period, the average density of the liq¬ 
uid deuterium in the target cell was thus determined to 
be p = (0.163 ± 0.002) g/cm 3 . 

The cylindrical portion of the liquid deuterium target 
cell measured 68 mm in diameter and 150 mm in length. 
The convex endcaps each added an additional 10 mm to 
the cell length at the central symmetry axis (which also 
corresponded to the photon-beam axis). Thus, the to¬ 
tal length of the liquid deuterium target cell along its 
symmetry axis was 170 mm. Due to the cell geometry 
and the divergent nature of the photon trajectories orig¬ 
inating from the radiator, the target thickness for each 
individual photon trajectory differed from this measured 
target length along the symmetry axis. A Monte Carlo 
simulation was thus employed to determine the effective 
target length for an “average” beam photon. The angular 
distribution of the photons emanating from the radiator 
was determined from a series of tagging-efficiency mea¬ 


surements taken with photon-beam collimators of differ¬ 
ent diameters. 

The effective target length was determined to be L = 
(166 ± 2) mm. re e ff was thus determined to be (8.10 ± 
0.20) x 10 23 (nuclei/cm 2 ), where the systematic uncer¬ 
tainty is the sum, in quadrature, of the uncertainties in 
p and L together with an additional 2% uncertainty to 
account for any potential misalignment of the symmetry 
axis of the target relative to the trajectory of the photon 
beam. 


C. Corrections 

1. Rate-dependent Corrections 

The rate-dependent correction arises from high elec¬ 
tron rates producing dead time in the single-hit FP TDCs 
and scalers. The nominal average electron rate of 1 MHz 
per FP channel, but the instantaneous rate fluctuated 
markedly and was as large as 4 MHz. The origins of these 
effects has been presented in detail in Ref. |36j, where a 
description of the Monte Carlo simulation of the FP elec¬ 
tronics used to quantify the effects is also presented. We 
note that this procedure has been used to successfully 
extract Compton scattering cross sections from carbon - 
see Refs. [ 3 !', ;38|. A short summary is presented below. 

The rate-dependent correction was the product of 
three terms 


/r — ,/ghost ./stolen/missed; ( 1 3 ) 

where /ghost was the correction arising from acciden¬ 
tal coincidences between the front and back scintillator 
planes in the FP hodoscope, /stolen corrected for prompt 
electrons that are not observed in the prompt peak be¬ 
cause an accidental electron stopped the single-hit TDC 
previously, and / m i sse d accounted for prompt electrons 
missed by the FP TDC but not the FP scaler due to 
deadtime in the electronics. Typical values for /ghost and 
/missed were 5% and 1% respectively, while /stolen ranged 
from 15 - 50% depending upon the FP rate for the runs in 
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TABLE III: Number of post-bremsstrahlung electrons, tagging efficiencies, and number of tagged photons. The relative sta¬ 
tistical and systematic uncertainties of N 7 are the same as those for e t ag- The upper four energy bins are from RP1 and the 
lower four from RP2. 


A 7 

(MeV) 

N e (x 10 13 ) 
[60°, 120°, 150°] 

Stag 

^^tag 

(statistical) 

^£tag 

(systematic) 

V 7 (x 10 12 ) 
[60°, 120°, 150°] 

69.6 

1.06, 1.06, 1.05 

0.421 

0.002 

0.005 

4.47, 4.47, 4.43 

77.9 

1.12, 1.12, 1.11 

0.423 

0.002 

0.005 

4.73, 4.73, 4.69 

86.1 

0.941, 0.941, 0.931 

0.426 

0.002 

0.005 

4.01, 4.01, 3.97 

93.4 

0.782, 0.782, 0.775 

0.425 

0.002 

0.005 

3.32, 3.32, 3.29 

85.8 

1.72, 1.72, 1.71 

0.456 

0.004 

0.007 

7.83, 7.83, 7.79 

94.8 

1.85, 1.85, 1.84 

0.459 

0.004 

0.008 

8.50, 8.50, 8.45 

103.8 

1.58, 1.58, 1.57 

0.460 

0.004 

0.008 

7.26, 7.26, 7.22 

112.1 

1.37, 1.37, 1.36 

0.458 

0.005 

0.008 

6.24, 6.24, 6.21 


question. A summary of the values of the rate-dependent 
correction /r is presented in Table IIVI where the first 
uncertainty is a scale systematic uncertainty common to 
all data points and the second is a kinematic-dependent 
systematic uncertainty that arises uniquely from / sto i en . 

2. Target-related Corrections 

The target-related correction was given by 

/t = jabs/fill/cell, (14) 

where / a b s was due to the absorption of beam photons by 
the liquid deuterium prior to scattering, fm was due to 
incomplete filling of the target cell, and / ce ii was due to 
beam photons scattering from the Kapton endcaps of the 
target cell. The absorption of beam photons by the liq¬ 
uid deuterium prior to scattering was determined using a 
Geant4 Monte Carlo which considered the effective target 
thickness discussed in Sect IIIIBI The correction / a b s was 
determined to be ~1.6%, and was known to better than 
3% relative uncertainty. 

During RP1, the liquid deuterium target did not fill 
completely. The liquid-deuterium level was observed 
each day by taking Polaroid images of the target cell 
which clearly showed the filled portion of the cell. Based 
on these images, it was determined that the top ~1.6 cm 
of the target cell was unfilled. In order to account for 
beam photons that passed through this unfilled portion 
of the target, a Monte Carlo simulation similar to the 
one used to determine K e ff was employed. The fraction 
of beam photons that struck the filled portion of the tar¬ 
get was determined, taking into account both the angular 
divergence of the beam photons and the uncertainty in 
the observed target-fill line. The correction fa u was de¬ 
termined to be 7%, and was known to better than 2.5% 
relative uncertainty. 

The contribution of the thin Kapton endcaps to the 
scattered-photon yield was investigated using a ~1 cm 
thick Kapton target. These Kapton data were subjected 
to the analysis detailed in Section IIII Al to extract the 
thick Kapton target yield Ikapton- The resulting correc¬ 
tion to the scattered-photon yield due to the thin Kapton 


endcaps was given by 

r I ^rkapton /, — \ 

J cell ) (I'-v 

where x was a factor used to scale the thick Kapton tar¬ 
get yield to the thin Kapton endcap yield. This scaling 
factor depended on the relative Kapton thicknesses, num¬ 
bers of incident photons, and target geometries which af¬ 
fected detector acceptances. The correction factors for 
DIANA (150°) and BUNI (120°) were consistent with 
unity within uncertainties. The average correction fac¬ 
tor or CATS (60°) was 94%. A systematic uncertainty 
of 2% in / C eii was determined for each detector. The 
target-related corrections are given in Table IIVI 

D. Uncertainties 

The dominant contribution to the statistical uncer¬ 
tainty came from the yield. The systematic uncertainties 
were divided into three classes: scale, angle-dependent, 
and kinematic-dependent. Scale uncertainties affected 
the data obtained at all angles and energies in a given 
run period equally. Angle-dependent uncertainties af¬ 
fected the results from the individual detectors differ¬ 
ently. Kinematic-dependent uncertainties affected each 
measured data point individually. We report these un¬ 
certainties separately. The sources of systematic uncer¬ 
tainties are listed in Table IV1 along with typical values. 

IV. RESULTS 
A. Cross Sections 

With the data analyzed as described above, we present 
the central result of our experiment in Table ED elastic 
Compton scattering cross sections on the deuteron. The 
results are also shown in Fig 0 along with those results 
at 66 MeV from Refs. |l6} and 94 MeV from Ref. [ 13 ] 
whose scattering angles are within 10° of ours. Statistical 
uncertainties only are shown on the data points. 
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TABLE IV: Rate- and target-dependent corrections. The first uncertainty in /r is a systematic that is common to all data 
points, while the second uncertainty is a kinematic-dependent systematic. The uncertainties in /t are systematic. The upper 
four energy bins are from RP1 and the lower four from RP2. 


E-y 

O 

O 

CD 


120° 


150° 


(MeV) 

/r 

/t 

/r 

/t 

/r 

/t 

69.6 

1.49 ± 0.05 ± 0.03 

1.03 ± 0.03 

1.60 ± 0.05 ± 0.02 

1.09 ± 0.04 

1.49 ± 0.05 ± 0.04 

1.09 ± 0.04 

77.9 

1.40 ± 0.04 ± 0.03 

1.03 ± 0.03 

1.49 ± 0.05 ± 0.02 

1.09 ± 0.04 

1.42 ± 0.04 ± 0.03 

1.09 ± 0.04 

86.1 

1.31 ± 0.04 ± 0.02 

1.02 ± 0.03 

1.36 ± 0.04 ± 0.01 

1.09 ± 0.04 

1.33 ± 0.04 ± 0.02 

1.09 ± 0.04 

93.4 

1.28 ± 0.04 ± 0.02 

1.02 ± 0.03 

1.32 ± 0.04 ± 0.01 

1.09 ± 0.04 

1.30 ± 0.04 ± 0.02 

1.09 ± 0.04 

85.8 

1.30 ± 0.04 ± 0.01 

0.96 ± 0.02 

1.58 ± 0.05 ± 0.02 

1.02 ± 0.02 

1.30 ± 0.04 ± 0.02 

1.02 ± 0.02 

94.8 

1.26 ± 0.04 ± 0.01 

0.96 ± 0.02 

1.48 ± 0.05 ± 0.01 

1.02 ± 0.02 

1.26 ± 0.04 ± 0.02 

1.02 ± 0.02 

103.8 

1.21 ± 0.04 ± 0.01 

0.96 ± 0.02 

1.39 ± 0.04 ± 0.01 

1.02 ± 0.02 

1.22 ± 0.04 ± 0.01 

1.02 ± 0.02 

112.1 

1.19 ± 0.04 ± 0.01 

0.96 ± 0.02 

1.34 ± 0.04 ± 0.01 

1.02 ± 0.02 

1.20 ± 0.04 ± 0.01 

1.02 ± 0.02 



E [MeV] 

FIG. 8: (Color online) Measurements of the deuteron 
Compton-scattering cross section from the current experiment 
(RPl [•] and RP2 [■]) are shown. Previous measurements 
from References [lf| (o), [il| (□), and El (o) are included 
for comparison. Statistical uncertainties only are shown. The 
solid and dashed lines represent the free and BSR-constrained 
fits to the present data only. The shaded region is obtained by 
varying the BSR-constrained fit for a s and /3 S by its statistical 
uncertainty. 


TABLE V: Sources and magnitudes of systematic uncertain¬ 
ties for the two run periods. 


Type 

Source 

RPl 

RP2 

Scale 

Tagging Efficiency 

1.1% 

1.6% 


Target Thickness 

2.5% 

2.5% 


Missed Trues 

1.5% 

1.5% 


Stolen Trues 

1.0% 

1.0% 


Ghost Events 

2.4% 

2.4% 


Kapton Cell 

2.0% 

2.0% 


Target Fill Level 

2.5% 

N/A 

Total 


5.2% 

4.7% 

Other uncertainties 


Statistical 

8-17% 

5-12% 

Angular 

Detector Acceptance 

3-4% 

3-4% 

Kinematic 

Yield Extraction 

2-11% 

2-11% 


Stolen Trues 

1-3% 

1-3% 


The current world data set of deuteron Compton 
scattering cross sections is comprised of three measure¬ 
ments |15l4l7| . The data reported here have uncer¬ 
tainties comparable to the previous data at low ener¬ 
gies (E < 70 MeV) and energy-bin widths considerably 
smaller than the previous high-energy (E = 95 MeV) 
measurement (Table IVTT1) . This experiment has doubled 
the number of data points in the world data set, in ad¬ 
dition to providing the first data above 100 MeV. As 
Fig. [5] indicates, all data sets are in excellent agreement 
within their respective statistical uncertainties. This is 
corroborated by a complementing angle transect at 66 
and 94.5 MeV, where these preceding measurements pro¬ 
vide additional data beyond our scattering angles; see 
Ref. El] f° r plots. Data consistency in the overlap is al¬ 
ready a strong indication that our results can well be 
embedded into the consistent world data set, extending 
it to higher energies. 

The figures also point to issues with two of our new 
data. The two points around 94.5 MeV and 60° are sep¬ 
arated by ^3<r, with the one from RP2 well in agreement 
with a SAL measurement and better agreeing with the 
overall trend. Another, more subtle, compatibility issue 
with the (112.1 MeV, 120°) point of the second run arises 
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TABLE VI: Measured cross sections for deuteron Compton scattering at the lab angles listed. The first uncertainty is statistical, 
the second is the overall normalization (5.2% first run period, 4.7% second), the third is due to detector acceptance, and the 
last is kinematic-dependent. The upper four energy bins are from RP1 and the lower four from RP2. 


E 7 |r(60°) &(120°) §7(150°) 

(MeV) (nb/sr) (nb/sr) (nb/sr) 


69.6 

15.7 

± 

2.6 

± 

o 

bo 

± 

0.7 

± 

1.0 

12.4 

± 

2.2 

± 

0.6 

± 

0.5 

± 

o 

bo 










77.9 

14.7 

± 

2.0 

± 

0.8 

± 

0.6 

± 

0.9 

15.0 

± 

1.3 

± 

0.8 

± 

0.6 

± 

0.3 

18.4 

± 

2.5 

± 

1.0 

± 

0.6 

± 

1.5 

86.1 

11.9 

± 

1.4 

± 

0.6 

± 

0.5 

± 

0.3 

15.7 

± 

1.4 

± 

0.8 

± 

0.7 

± 

0.3 

15.7 

± 

2.3 

± 

0.8 

± 

0.5 

± 

0.8 

93.4 

8.1 

± 

1.2 

± 

0.4 

± 

0.3 

± 

0.3 

16.0 

± 

1.3 

± 

0.8 

± 

0.7 

± 

0.4 

13.7 

± 

2.2 

± 

0.7 

± 

0.4 

± 

1.5 

85.8 

13.8 

± 

1.7 

± 

0.6 

± 

0.6 

± 

1.5 

13.4 

± 

1.0 

± 

0.6 

± 

0.6 

± 

1.2 

16.8 

± 

2.0 

± 

0.8 

± 

0.5 

± 

0.7 

94.8 

15.4 

± 

1.5 

± 

0.7 

± 

0.7 

± 

1.2 

14.1 

± 

OO 

o 

± 

0.7 

± 

0.6 

± 

0.5 

15.1 

± 

1.7 

± 

0.7 

± 

0.5 

± 

0.5 

103.8 

11.9 

± 

1.1 

± 

0.6 

± 

0.5 

± 

0.3 

11.8 

± 

0.7 

± 

0.6 

± 

0.5 

± 

0.4 

15.7 

± 

1.6 

± 

0.7 

± 

0.5 

± 

0.8 

112.1 

8.8 

± 

1.0 

± 

0.4 

± 

0.4 

± 

0.2 

9.8 

± 

0.7 

± 

0.5 

± 

0.4 

± 

0.2 

13.0 

± 

1.5 

± 

0.6 

± 

0.4 

± 

0.4 


TABLE VII: Comparison of reported results from measure¬ 
ments of deuteron Compton scattering. 


Ref. 

E 

(MeV) 

AE 

(MeV) 

Number 
of Points 

Normalization 

Uncertainty 

[15| 

49,69 

6.5,7.7 

6 

4% 

[III 

95 

21 

5 

5-6% 

m 

55,67 

10 

18 

6-14% 

This work 

70-112 

7.3-9.0 

23 

5% 


when comparing the data to fits, see also Refs. (l. 23j. 
For both instances, the discrepancies could not be traced 
back to specific data or analysis issues. We therefore 
report these points but flag them as potential outliers, 
presenting results with them included in our extraction 
of the isoscalar nucleon polarizabilities in this paper. 

B. Nucleon Polarizabilities 

As argued in the Introduction, the goal of this exper¬ 
iment was to provide new, high-quality data in support 
of an extraction of the neutron polarizabilities. Such a 
determination of these two-photon response parameters 
needs theoretical input, but it will also allow for a better 
understanding of the quality of our data. 

First, the effects of nuclear binding and charged meson- 
exchange currents inside the deuteron must be sub¬ 
tracted. They contribute a significant fraction of the 
deuteron cross section at the energies we measured [39| . 
as well as in the zero-energy limit in which aEi and 
13mi are defined [40] . For the individual amplitudes, one 
then needs to subtract the Powell amplitudes of point-like 
spin -!/2 particles with anomalous magnetic moments. 

More importantly, however, the nonzero-energy data 
must be related to the static point. In a simplistic ex¬ 
trapolation, one may be tempted to use a cross-section fit 
which identifies the electric and magnetic polarizabilities 
as angle-dependent contributions to the terms quadratic 
in the photon energy. This is, however, not permissible 
since all of our new data, and most of the world data, 
are well beyond its realm of applicability. At energies 


at and above 100 MeV, the energy-dependent effects of 
the pion cloud and of the A(1232) excitation are impor¬ 
tant, and the pion-production threshold induces the first 
non-analyticity in the single-nucleon amplitudes. A low- 
energy expansion proceeds thus in powers of ui/m„ and 
becomes quickly useless. 

A viable low-energy parametrization should hence con¬ 
sistently account for all these effects and provide a sys¬ 
tematically improvable estimate of residual uncertainties. 
Chiral Effective Field Theory (%EFT) is ideally suited 
for this task. It model-independently encodes the correct 
symmetries and effective degrees of freedom of low-energy 
QCD, with a small dimensionless parameter to systemat¬ 
ically improve the description of higher-order effects. It 
is well established that %EFT predicts the energy depen¬ 
dence of the single-nucleon Compton scattering response 
over the full range of data with high accuracy, includ¬ 
ing spin effects Q, and that it consistently accounts for 
nuclear binding as well as meson-exchange and nucleon- 
nucleon rescattering effects. All these aspects are neces¬ 
sary to restore the Thomson limit on the deuteron. 

We therefore turn to the yEFT description which was 
used in previous high-accuracy descriptions of proton and 
deuteron Compton scattering. As its ingredients at next- 
to-leading order in ckei and (3mi (order e 2 S 3 ) have re¬ 
cently been described summarily, we refer to Sect. 5.3 of 
Ref. [2| for details. The interactions between nucleons, 
pions and the A(1232) resonance arc fully determined ex¬ 
cept for the two which parametrize short-distance contri¬ 
butions to the scalar polarizabilities. The dependence of 
the resulting extraction of a and /3 on the deuteron wave 
function and NN potential was shown to be negligible, 
and residual theoretical uncertainties were estimated as 
±0.8 canonical units. In the results reported here, we 
use the same fit procedure and parameters as reported in 
Ref. As described there, we add kinematic-dependent 
and angle-dependent systematic uncertainties in quadra¬ 
ture to the statistical uncertainty, and treat correlated 
overall systematic uncertainties by a floating normaliza¬ 
tion. This determination is based on the entire data set 
of our experiment, but does not include the other world 
data. Finally, we treat the two run periods as statistically 
independent data sets. Treating them as a single data set 














12 


TABLE VIII: Isoscalar nucleon polarizabilities and \ 2 P er 
degree of freedom extracted from the world elastic deuteron 
Compton scattering data set prior to our data, from our data, 
and from the new world data set including our data, each 
using both a “free” fit to a s and and a determination 
constrained by the BSR. Uncertainties are statistical only and 
anticorrelated for the BSR-constrained fits. 


set 

constraint 


Ps 

X 2 /d.o.f 

Ref. 

old world 

free 

10.5 ±2.0 

3.6 ±1.0 

24.3/24 


BSR 

10.9 ±0.9 

3.6 =F 0.9 

24.3/25 

—ra- 

this work 

free 

13.2 ± 1.4 

3.2 ± 1.1 

40.6/19 


BSR 

12.1 ±0.8 

2.4 =f 0.8 

41.7/20 


new world 

free 

BSR 

11.1 ±0.9 
11.1 ±0.6 

3.3 ±0.6 

3.4 =f 0.6 

49.2/43 

45.2/44 



changes the following conclusions at most marginally. 

Table IVIIII summarizes the findings of our fit, in the 
context of the previous extraction of the polarizabilities 
and the determination based on the new world data set 
which includes our data but with the two points previ¬ 
ously mentioned discarded - these results being already 
reported in Ref. jlj and Eq. ©■ For our extraction, as 
for the others, the values of the independent fit to a s and 
Ps demonstrate excellent consistency of a s + p s with the 
isoscalar BSR, Eq. ©• The BSR can therefore be used 
to reduce the number of parameters from two to one, de¬ 
creasing the statistical uncertainties. In the fit to only 
the new data, the normalization of RP1 floats by about 
2% down, and that of RP2 by 3.5% up, i.e. well within 
the overall correlated uncertainty. The statistical uncer¬ 
tainties are smaller than for the previous world data set, 
but the x 2 per degree of freedom is larger. As hinted 
above, this can be attributed to two points which be¬ 
tween themselves contribute about 20 units to x 2 > while 
changing the central values only within the statistical un¬ 
certainties. A complementary publication [23j provides 
details in the context of the construction of a consistent 
database. Here, we reiterate that a careful analysis of our 
data-taking and analysis procedures showed no intrinsic 
experimental reason why these points should be special. 

Our central values for either fit agree with those of the 
old and new world database extractions well within the 
systematic uncertainties only, not accounting for theoret¬ 
ical and BSR uncertainties. Since it effectively doubles 
the world data, it is therefore no surprise that the new 
world average reported in Ref. [l ( ] is hardly shifted but 
its statistical uncertainty is reduced to l/-\/2 ss 70% of 
the previous one. 

With only our current data, we thus obtain the final 
values of the BSR-constrained fit as 

a s = 12.1 ± 0.8(stat) ± 0.2(BSR) ± 0.8(th) 

P s = 2.4 =f 0.8(stat) ± 0.2(BSR) ± 0.8(th), 

where the uncertainties from theory and from the BSR 
constraint are listed separately. One finally extracts the 


neutron polarizabilities by combining with the proton 
values of Eq. © which were determined using the same 
xEFT approach and fit philosophy 

a n = 13.55 ± 1.6(stat) ± 0.2(BSR) ± 0.8(th) 
p n = 1.65 =F 1.6(stat) ± 0.2(BSR) ± 0.8(th). 

We list these values for completeness only, since those 
obtained from the new, statistically consistent world data 
set p]] supersede them in accuracy 

a n = 11.55 ± 1.25(stat) ± 0.2(BSR) ± 0.8(th) 

P n = 3.65 =F 1.25(stat) ± 0.2(BSR) ± 0.8(th). 

In conclusion, careful statistical tests indicate that 
our measurement provides new, high-quality data whose 
analysis is well-understood. 


V. SUMMARY 

This paper reports new Compton-scattering cross sec¬ 
tions for deuterium over an energy range from 70-112 
MeV, at angles of 60°, 120° and 150°. The data points 
are in excellent agreement with previously published re¬ 
sults. An analysis using %EFT extracts isoscalar nucleon 
polarizabilities that agree, within uncertainties, with pre¬ 
vious extractions. This measurement represents the first 
new result from deuteron Compton scattering in more 
than ten years, nearly doubles the number of data points 
in the global data set, and extends the maximum energy 
by almost 20 MeV. Furthermore, this data set reduces 
the statistical uncertainty of the extracted values of a n 
and p n by a factor of ~l/3. 
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